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. Abstract 

. In this paper, we study a linear system related to the 2d system of Eu- 

ler equations with thermal conduction in the quasi-isobaric approximation 
of KuU-Anisimov [14]. This model is used for the study of the ablation 
front instability, which appears in the problem of inertial confinement fu- 
sion. The heat flux Q is given by the Fourier law T~'^Q proportional to 
VT, where v > 1 is the thermal conduction index, and the external force 
^ ' is a gravity field g = —gCx- This physical system contains a mixing region, 

' in which the density of the gaz varies quickly, and one denotes by Lo an 

0^ ■ associated characteristic length. The fluid velocity in the denser region is 

' denoted by Va- 

O 

. The system of equations is linearized around a stationary solution, and 

' each perturbed quantity u is written using the normal modes method 

!>■ ; 

O ■ u{x, z, t) = 5R(u(a;, fc, ^)e»'=-+^v^*) 

in order to take into account an increasing solution in time. 
, The resulting linear system is a non self-adjoint fifth order system. Its 

» I ■ coefficients depend on x and on physical parameters a, /3, a and f3 being 

' two dimensionless physical constants, given by q/3 — kLo and ^ = 

(introduced in [5]). We study the existence of bounded solutions of this 
system in the limit a 0, under the condition /3 £ [/3o, and the 
assumption 5R7 G [Oi 111 — ]^ (regime that we studied for a simpler 
model in [5]) calculating the Evans function Ev{a, f3,^) associated with 
this linear system. 

Using rigorous constructions of decreasing at ±00 solutions of systems of 
ODE, we prove that, for /? G [/?o, 3^7 G [0, ^], I7I < there exists 
ai > such that there is no bounded solution of the linearized system for 

< Q < Ql. 
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Necessarily, for any M > and /3o > there exists ai > such that, for 
< a < ai and /3 £ [/3o, ^], an admissible value 7(0:, /3) such that there 
exists a bounded solution of the linearized system satisfying I7I < M is 
such that 5R7 ^ [0,M]. 
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Introduction 



This paper is devoted to the precise calculus of the Evans function Ev{a, P^-f) 
of the normal mode formulation of the linearized system of equations associated 
with the quasi-isobaric (low Mach number) model. The calculus of this Evans 
function is not classical, because the matrix of the differential system has sin- 
gular coefficients, and because these coefficients do not behave exponentially in 
the spatial variable. However, usual techniques of ordinary differential equa- 
tions and introduction of a Fuchsian problem allow us to calculate this Evans 
function under certain assumptions on the parameters a, (3, 7 introduced in the 
Abstract. 

In this Introduction, we first describe the physical model ([5]), define what is 
called a linear growth rate of the linearized system associated with a stationary 
solution of this physical model, then finally describe the contents of each step 
of the proof of the main Theorem (Theorem 2]) . 



0.1 Physical model 

We consider a compressible fluid characterized by its density p, its velocity (u, v) 
and its temperature T in a gravity field g = — \g\ex- We assume that this fluid 
has the following properties: 

a) when x goes to +00, for all z we have p — > pa, iu,v) — > (—14,0) and 

b) the functions {p, u, v,p, T), where p is linked to the pressure in the fluid, 
satisfy the system of the Euler equations in two dimensions (x, z) with thermal 
conduction in the quasi-isobaric approximation for a perfect gaz: 

div(Cppru-|-(3) =0, (1) 

where Cp is the calorific capacity of the fluid, the heat conduuction flux Q being 
given by the Fourier law 

Q - -k{T)VT, (2) 

the thermal conduction law is 

k{T) = i^oT^ (3) 
where v is the thermal conduction indice. 

We introduce a characteristic length Lq associated with the thermal properties 
of the fluid 

" = r nV ^' 

Physical values of Lo for the case of the ICF arc of order 10^*^ meters. 

Under a quasi-isobaric assumption, the system modelizing the ablation model 
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was given by H.J. Kull [13] and appears for example in P. L. Lions [T7]. It states 

(5) 



dtp + div(pu) = 
dt{pu) + div(pu ®u + p)=pg 

pT = PaTa 

Q = -k{T)VT 
[ div(CppTii + Q) =0. 



This model can be derived either from the low Mach approximation (Majda [T^], 
Dellacherie ^) or the quasi-isobaric approximation (Kull [14^, KuU-Anisimov 
[15], Masse [20]). See a short analysis in Section[TJ A stationary laminar solution 
of the system© is {p[){x),uq{x),Q,P[){x),Tq{x)), where we introduce a function 
^ such that 



Pa{x) = Pai{ — ),uo{x) 

1^0 



PaVg ^ Vg 

'po{x) a-t) 



andpo(a;) satisfies 



Poix) 



P^V^ 

ra ^ a 

Po{x) 



po{s)ds = poixoy I , ^ 

x„ PO(X) 



_ PaTg _ Tg 

PlVg' 



The function ^ is the solution of the differential equation 

ay 

such that ^(0) = Note that, in this case 



(6) 



Po{Loy) =Po(0) 
Introduce 



I'+l 



PaVg - PaV^ 



9L0 

Va' JO 



C{t)dt] 



Zip) 



{v + l)p''+^ 
The system of unknowns that we consider is 

/ pu 

U = 



(7) 



pw +p 
puv 
Z{p) 

\ u-LoVad^{Z{p)) J 



From U, we recover p from Z{p), u ~ v — and p = p -\- pu^ ^ . 

For this choice of unknowns, we introduce Fi{U) — U. There exists three 
explicit functions Fq, F2 and F3 such that the system (O is equivalent to the 
system on U 



dtFo{U)+d,Fi{U)+d,iF2iU,d,U)) = FsiU) 

z-HUi) 



(8) 
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A stationnary laminar solution Uq of this system depending only on x is Uo{x) 
such that 

±(F,iUo{x))) ^ FsiUoix)). 

The identity FiiU) — U\s the natural choice when one studies a basic solution 
depending only on the variable x. 

0.2 Definition of a linear growth rate 

We linearize ([8|) around Uo{x). Denote by U the unknowns U — Uo{x). The 
linearized system writes (Vi and V2 denotes the gradient of F2 with respect to 
the first and second set of variables U and dzU): 

d,U+VFo{Uo{x))dtU+ViF2iUo{x),0)d,U+\72F2{Uo{x),0)dl2U ^\7F3{Uo{x))U 

(9) 

which can be rewritten M{x,dx,dz,d'^2,dt)U — 0. Note that its coefficients 
depend on x through the stationnary solution. 

We are now ready to introduce the definition of a linear growth rate for a 
non linear system around a laminar solution: 

Definition 1 Let M{x,dx,dz,d^2,dt)U = be the linearized system 

We call a linear growth rate of this system for the wave number k a value 
of a ( depending on k) such that 5Rct > and there exists a non-trivial solution 
U{x, fc, cr) of the system 

M{x,-^,ik,-k^,c7)U{x,k,a) =0 (10) 
such that U is bounded and going to when x goes to ±00. The function 

U{x,k,a)e'^''e^* 
is called a normal mode solution of the system. 
The normal mode system associated with © is 

'^+VFo{U^{x))aV + ikViF2{Uo{x),Q)V -k-'V2F2{Uo.m ^^F^iM^W- 
ax 

(11) 

The scope of this paper is to find bounded non trivial solutions of (|lip . and 
associated values of a if any. If such a solution exist, it will lead to a normal mode 
solution of the linearized system. Note that, in the set-up we described, different 
physical parameters appear, namely fc, Lq, Va, g- As the classical growth rate of 
Rayleigh is equal to ( ^^~^^ fffc)2 for the discontinuity model [22], and as we 
proved ([5], [12]) that this value was the limit of the growth rate when kL^ goes 
to zero, we are led to introducing the following quantities 

e = A:Lo, Fr = ^ (12) 
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and 

(13) 

The aim of this paper is to study the existence of a growth rate 7, 5R7 > 0, in 
the hmit Ln ^ Q when the Froude number Fr is of order , which means that 

a 0,/3 > 0. 



Remark Other regimes rely on different assumptions on a and we refer to 
[13] for the resuhs that can be obtained for this model in the high frequency 
regime k +00. In this other regime the scaling writes 

£ large, e^Fr < C , where C is a constant. 
In the first section, we study the physical origin of the model and derive the 
properties of the stationary solution, where the associated density profile satis- 
fies: 

Pa{x) —>■ when x — > —00 
pa{x) -pa — Ce^^ ,x +00 
po{x)\x\^ Pa{^)",x-^ -00. 

We then derive the linearized system, which is a fifth order differential system 
whose coefficients depend on po{x) and are singular when poix) go to zero. Note 
that this is not a classical case for the study of such systems and that this leads 
to rather tricky methods. 

In the second section, we recall the general set-up for the calculation of the 
Evans function of the linearized system, and we give the induced differential 
systems in A^^{K^) for n — 2,3. Note that the field is IR for real values of 
7 and K ^ C for complex values of 7. This Evans function is (related to) the 
vectorial product of the normalized solution in A^{K^) which has the greatest 
decay when x — > +00 and of the normalized solution in A^{K^) which has the 
greatest decay at a; ^ — cx). 

In the third section, we identify the solutions of the system in A'^{K^) deduced 
from (jlOp for x —^ +00. In this region, we use the exponential behavior of the 
profile to obtain the classical analytic expansion of the normalized solution of 
the system in A^(if^). There exists ^0 G]0, 1[ (corresponding to yo G R through 
^iuo) = Co) such that this analytic expansion is valid for C(x^) ^ Coi that is 
X > LqUo. Note that, however, the expansion of this solution cannot be obtained 
by the techniques developed in Zumbrun aZ [4], because the Gap lemma as- 
sumptions are not fulfilled. 

A general feature in the calculation of the Evans function is to obtain an over- 
lapping region of definition between the solution well behaved at -f 00 and the 
solution well behaved at —00. A first step to achieve this overlap is then to 
prove that there exists ao > and R > such that, for all a < ^ < the 
solution obtained for x > LqJ/o can be extended in [^".^(q:, (^), Lqi/q] where 

X,(a,C) ^ i 1 
^0 
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This is the aim of the fifth section. The behavior of the solution in the region 
[ct~ Co " , Co] when a — > is different from the classical analytic expansion in a 
for ?/ e [?/0i +00 [ and it is the aim of Sections 13.31 and 13.41 

Once this extension is done, an easy calculus is the calculus of a growth rate 
associated with the following stationnary solution, characterized by its density 
profile, for a (0 such that (0 < 

P*(a;) = < i"i (14) 
\ Paa-Qo " . 

This calculus is an improvement of the discontinuity model of Piriz, Sanz and 
Ibanez and it is the aim of Section HI 

When the profile is not constant in the region ] — 00, X*(a, (^0)] (that is for the 
full model), the system leads to a fuchsian problem in the region x ^ —00, and 
we use the hypergeometric equation (see [16]). The solution of the system in 
tc'iK^^ deduced from pO|) is identified in any region of the form ^~'^-<~^^ 
for every ioj which means that a; s] — 00, -~'^\- The results of the analysis of 
these solutions is summarized in Theorem [1] 

The study of the roots of the Evans function is the aim of Section [5] and we 
summarize the method here. From the relation (,(y)\y\~ ^ v~~ when y —00, 
we deduce that —aPX^{a) > when a ^ 0. Hence for all < Oa < ao 

there exists to > such that the regions ] — 00,—^] and [ ^^IT"* ' ^0]' '^^ere 
C(yo) = Co, overlap. 

We then express the Evans function Ev{a, P^-f) of the system at a point of 
["o^'^T^]' "^^^ limit when a ^ and to small exists and we write its 
expression in terms of r = ^, /3 and to > 0. As it does not depend on to we 
study the limit when to ^ +00, hence proving that the only positive value of r 
which is admissible is r = 1. We deduce a contradiction, proving that there is 
no growth rate ( of positive real part) for the system. This can be stated as 



Theorem Let M be given. There exists a* > such that, for < 

a < a*, /3 e [jg, AI], the Evans function Ev{a, /?, 7) of the system has no 
root for I7I < M, 3*7 e [0,A/]. 



1 Derivation of the quasi-isobaric model 
1.1 The physical approximations 

The general equations are the thermal hydrodynaniic equations, written in a 
non conservative form: 

dtp + div (pit) — 

dt (pu) + div(pu <S)U+p) = pg (15) 
p{dt + u.V)h - {dt + u.V)p ^ -div(g + Ji). 
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where Cp and C„ are the classical tehrmodynamic calorific capacities at constant 
pressure and at constant volume, h is the enthalpy h = CpT, the pressure 
and the density being given by the equation of state p = {Cp — Cv)pT, Q = 
—k{T)WT, Ji = (in our assumption the energy given to the system is 0). The 
quasi-isobaric approximation writes 

Cp - Cy Sp ST 
— « — . 

Following L. Masse [20], this relies on the two hypotheses 

M^«l,— «1. (16) 

Hence the quasi-isobaric model relies on a low Mach hypothesis. 
The equation of the energy rewrites 

aivu H -divQ = . 

ph ph 

This equation is approximated by ([T]) as we will see below. What follows is a 
formal derivation of the quasi-isobaric model under a low mach hypothesis. It 
is closely related to the method given by Majda ,191, Dellacherie [B]. 



1.2 Adimensionnalization and low Mach expansion 

Use the reference density pa, the reference velocity Va, and the reference pressure 
Ps associated with the sound velocity Cg such that c^^Pa = Ps- The Mach 
number is thus M — —. 

Write p = PaP' , u = Vail' , p — PsP' ■ The system of equations (jl5|) rewrites 

f V-^dtp' + div{p'u') ^0 

V-^dt{p'u') + div(p'^' ®u' + = 

If we assume that all the quantities have an asymptotic expansion in M, in 
particular p' — p'q{x, z, t) + M^p{x, z, t, M) we have the following relations from 
the momentum equations 

dxPa = 0, dzPo = 

hence p'q depends only on t. This is the same result as in the analysis of Del- 
lacherie [B]. In this model, we assume that the pressure p'q is constant, because 
we assume that the ground state for the equations is stationnary. 

Replacing the relation p'{x, z, t, M) = p'q + M^p in the energy equation we 
obtain 

" Cpip'o + M^p) ^ p'o + M^p 
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Finally, using Q = Kq ( (^^itcX' ^""^ (S-cC)p' ' deduce that 

hence the formal analysis leads to the equation 

div{V^'u' + LoVZ{p)) = 0{M^) 

where we used ^ = deduced from the relation p = [Cp — Cy)pT. The 
resulting equation can be written div(CppTu + Q) = 0, pT = paTa hence (fT|). 

Finally, in the momentum equations, rewriting — + P ^-^d using Pq 
constant, we obtain the equations 

V-^dt{p'u') + div(p'w' ® u'+pid) = -^p'. 

Note finally that the relation ([1]) and the relation pT — paTa lead to the equation 
on p: 

{dt + u.V)Z{p) -{!, + l)LoVaZ{p)AZ{p) - 0. (17) 

1.3 Study of the stationnary solution 

The resulting system of equations that models our phenomenon is thus 

dtp + drv{pu) = 

dt (pu) + div{pu (g) M + pid) ~ pg (18) 
div{u + LoVaVZ{p)) = 0. 

A stationnary laminar solution satisfies 

Po{x)uo{x) = -paVa 
_d_ 

If 



£{Poix)uo{x)'^ +Po{x)) = ~-po{x)g 



dx 



{uo{x)-LoVaZ{po{x)))^0. 



where the first relation is a consequence of the mass conservation equation, the 
constant pqUq being identified through its limit at a; — *■ +oo. Hence uo{x) = 



PaVa 



leading to the equation on po: 



^LoVa— Pq{x) p„+ ^Cq. 



Poix) dx 
As y = this equation rewrites 

dy Va^ ■ 

If Co = 0, the equation becomes ") — —^i hence ^"'^ = — vy hence ^ 
is not defined for y > —. We cannot consider this solution. 
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If Co > 0, ^ is increasing, hence if it is majorated, it has a Umit I > when 
y +00, this limit I satisfies l'^'^^ + ^l'^^^ = 0, which is impossible. We deduce 
that Co < 0, hence the equation is 

dy Va 

hence from the resolution of the equation we deduce ^ l^T' as ^ — > 1, 
|Cq| — Va and the resulting equation is jH]). 

This equation has a unique constant solution ^ = 1. The low Mach approxi- 
mation of S. Dellacherie |6j for a bubble model uses this stationnary solution as 
base solution. When we consider a non constant solution, we have the following 
Lemma which gives the behavior of the solution when x goes to — cxd. Introduce 
the function h^^j such that = Jq :;p^TiT(iZ;^ 1 a-i^d n — [v], = ^* 

and y, = Y^^^^ ^C'" ~ 

Lemma 1 There exists > such that for < t < t^ there exists a unique 
continuous function g{t), g{0) — 1 solution of 



mr = ~t''9itnh4tg{t))+yo)+Y.t^igit)) 



p " 



v — p 

p=0 ^ 
The function g has a Taylor expansion at t = 0. 

1. Behavior when y — > —00 

For y < —jYy^ J ™6 have the identity 

There exists a function r{t, e) such that 

-ia--)r = -+eh-'-ir{t,e). (19) 
e e vt 

There exists > and eo > such that r{t, e) is bounded for t > tQ,0 < 
£ <£o 

2. Auxiliary function at —00: We introduce S{t, e) — — J^°° s^^^ir{s, e)ds. 
This function satisfies S{t, e)ti uniformy bounded for t > to and < £ < 
£o- 

3. Behavior when y — )■ +00 

We have 1 - £_{y) = C{y)e-y, with C{y) ~* exp{yo - i - .. - -^j^ + 

fl (l-i)'"~")drj N 
Jo r/''-"(l-r;) ^■ 
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The proof of this Lemma uses the following relation 

n 

* — n ^ 



p=0 

The equality yields 



1 

^ p=0 

Introduce S, = tg and t = (— ^)^, such that 2/ = — We obtain the equality 



Wc introduce ^{t,g) = g'' - iT^t^g" - ^t" g''[h{u}{fg) + yo]- Wc have 

\l/(0, 1) = and 9g\l/(0, 1) = v > Q hence by the implicit function theorem, 
there exists a unique solution of ^{t,g) = in the neighborhood of (1,0). 

The behavior of the solution is more classical in the neighborhood of +oo, 
thanks to the equality 

In 1 - ,e = yo - V 7 — — + / ^ - y 

^ ^ Jo »7 "(1-??) 

hence 

with C(y) ^ yo - Ep=o + Jo^ ^^^"^ ^ ^ "'■'^^ Lemma is 

proven. 

The particular case u = 2.5 In this case, explicit calculations lead to the 
following implicit relation 

Hence, for y = — | and ^ = rjei we get 

2 2 4 ^ , 1 + eK/ri 
t = — ^{l+s^rj + s'sr]^) - eC- eln 

2 1 — e 5 

We write T]{e,t) = e). Considering the limit, for t > fixed, in the 

previous equality we obtain 

\ime^og{t,e) = 1 
and we construct step by step the expansion of g{t,e) in e. 
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1.4 Physical interpretation of the model 

We insist finally on the fact that this system of equations is only a theoretical 
model: the equation satisfied by po is 

hence the pressure is not bounded. The total pressure writes 

P{x) = [Cp - OpaTa + M^poiLox) 

which leading order term is constant and low order term in Mach is not bounded. 
Hence a good way of calling this model could be to call it a relative isobaric 
model or a relative low mach model. See Majda and Sethian |19j, Embid 
[7], or P.L. Lions [TJ for other remarks on this modelling. 

Introduce the function M{y) such that M^{y) — ■'^r^: where Cs{y) is the 

sound velocity at a point of the fluid given by Cs{yY — §^f(fy' ^^^^ ^^^^ 
call the Mach number. We have 

Lemma 2 The Mach number of the stationnary solution is bounded for v >2 
and y < —C. 

Proof As M^di^ - P-^.- - - Pa^g ^ 

i-roor AS lu \y) - ^^(^y)c'iiy) " ip(y)i(y) ^ (Cp-C„)p„T„^(j,)+M^po(y)«y)M^ 

which has a finite limit when y — s- — oo under the condition ly > 2. More- 
over, as poiy) — oo when y +00, there exists a point where the pressure 
vanishes. There exists a constant C such that P{y) is bounded below > on 
] — 00, — C] and — C is a positive constant of order (3^) • 

The low Mach number assumption is relevant (in particular when looking at 
the temperature of ablation and the density in the ablated fluid). 



1.5 Linearisation of the equations 

To simplify the notations of what follows, we denote by / a quantity appearing 
in the Euler system of equations, by /o its stationnary leading order term, and 
by / the perturbation of order 1 normalized by the physical quantity pa for the 
density, paVa for the impulsion, and PaV^ for a pressure term. In the normal 
mode study in the vicinity of a proflle depending on x, it is pertinent to linearize 
the variables on which acts the derivative that is pu, p{u)^ +p and pud, 
which writes 

pu = -paVa + PaVaXi 

p{uf +p^ poix) + po{x)iuoix))^ + PaV^X2 (20) 

puv = -ipaV^xs 

The coupling with the equations of the energy is made through the pertur- 
bation of density p. Introduce 

f = ^ - LaVaV{Z{po)) - LoVa\7{Z{p) - Z{po)). 
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The thermal perturbed quantities are 



X4, = Z{po) - Z{p),X5 = — (n + Va). 



(21) 



With the choice of unknowns (xi, 2:2, 0:3, a;4, xs) (which corresponds to the un- 
knowns pu, pv? -\-p, puv, —Z[p), Ti), denoting by the inverse function of Z, 
Z~^{f) = Pa{{v + 1)/)'^, the non linear system (|18p is equivalent to: 

' V-^dt{Z-\Z{po) - x^)) + pad^xi + id,{Z-\Z{po) - xi)j^) = 
V~^dtxi + PadxX2 - ipadzXz = -^{Z^'^{Z{po) - 2:4)) + ^Po{x) 

V-^dt{Z-HZ{p^)-X4)^)-pakx^ 

-I-/, r) (t^ -I- Z ^(Z(pQ)~Xi) Pa Pa{-xif N _ p. 

+PaC'4.T2 + + ^ - Z-^(Z(po)-x,)l - ^ 



Pa(l-Xl) 



LodxX4+ z-i(Z{po)-Xi) po(x) 

dxX5 + dz{LodzXi + j = 0. 



2^5 = 



(22) 

In the system (|18p . one needs to obtain the linearization of pv and of p + p{vY 
in terms of xi,X2,x^. It is a consequence of 

« = ,P = Pa 4 + P ,P ^X2 — h -,w = -— — - 

hence the approximations (dropping the terms of order 2 at least) 

V = iVaX3,p Po + PaV^{x2 + + T^),"" = WO + "^P- 

Moreover, the identity a;4 = Z{paS,) — Z{pa{£, + pj) leads to x^ — — Z'(po)PaP + 
0{p'^), hence 



^4 = ^+0(^). 



Moreover, as 



Tl It 



- K - L„Vadx{Z{pa{^ + p)) - Z{paO) 



one deduces 

2^5 = Kr^(" - Mo) - PaLt)dx{Z'{paOp)- 

The linearized system is 



f F-'9tp + a.xi+e9.(«a.3) =0 
V~^dtxi + dxX2 - dziixs) = ^p 
V-^dt{i^X3) ~ 9,(«a;3) + r'd.\2xi - TV) = 
p = -(pa^'(PaC))"^a;4 

LodxiPaZ'iPaOp) +X5 + |f P = 

, 9^X5 + dz{LodzX4 + ixa) = 



(23) 
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Using the new variable y such that x = Loy, one deduces the relation 



X5 dyX/\. 



p + ^{y)xi „ , cu , xi 

~ OyXA + E, Xi ' 



' i{ymy) + P) 

Write the following normal mode expression: 



o{ 



xip + 



( ^1 \ 

iX3 
\ X5 J 



( 3i, \ 



^2 



(24) 



V 25 / 



Assume 



/ £1 \ 



^ikz+-/V9^t is a solution of the linearized system. This rewrites 



is solu- 



V 3l5 / 

as a system of ordinary differential equations on 



tion of this linear system of ODE, then {xi,X2,ix3,X4,X5y is solution of ([23 
Note that, in this case ix^ is real. 



21 






/ 


2l 


\ 










£2 








. If 








■£4 








■£4 




^5 






\ 


^5 


/ 



Remark on complex growth rates Note that, if 7 is complex, the solution 

X.2 

ix^ is also complex and depend on y, k, 7. More precisely, introduce yj, Zj 
such that Xj = yj + izj. We have thus 



5R( 



/ 3Ii \ 

X.2 

V £5 / 



( j/i COS fcz — zi sin fcz ^ 
2/2 COS kz — Z2 sin kz 
— Z3 cos kz — 2/3 sin kz 
2/4 cos kz — Z4 sin kz 
y 2/5 cosfcz — Z5 sin/cz j 



and 2/j and Zj are known through the decomposition of the initial perturbation 
in the even part and the odd part in z. 



Hence, from a complex solution [x_-^ 



; ■ • ■ I ;k5 



)*(a:, fc, 7) of the normal mode system. 



one deduces a solution of the perturbation system with a known initial condition. 
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Normal modes system The equations for the normal modes associated with 
the mass conservation and the momentum equation are 



+ if + = 



dx 



- a(3x3 + a^xi + |p = 
dm _ aj^xs + aP{x2 + fxi + ^) = 0. 



(25) 



The normal mode formulation of the linearized energy equation is 

— h ial3{—ix'i + iajSxA ~ 0. 

dy 

Let X be given by 



X 



( ^1 \ 

X2 

X4 
\ X5 J 



The linearized system on X is 

dX 



— + Mo{^,a,P,j)X^O 
dy 



(26) 



where the matrix Mq is given by 

/ 



2a/3 



V 






aP 





a/3e 


a/3 



~a^P^ 



\ 




-1 
/ 



(27) 



Prom now on, we will call this system the Kull-Anisimov system. The 

eigenvalues of —Mq are given by the classical resul10 (see [14 , 21 ), and help 
us to study the solution at ±oo: 



Proposition 1 • The eigenvalues o/— Mq are 

Ao(0 = a'y£„\a,+ = a/3, Aa,- = -a/3, A+(^), A_(^) 
where (the square root is chosen of positive real part) 



A±(0 



± 



v+l 



a2/32 



(28) 



(29) 



The eigenvalues Aq, Aa.± are called the hydrodynamic modes, the eigenval- 
ues A±(^) are called the thermal modes. 



change of unknowns (that is a general X = R(y)X) lead to a different set of eigenvalues, 
however if the matrix R{y) depend only on and is a function of § G [0, 1] then the 
limit of the eigenvalues when y — > ±oo is the same as the limit of the eigenvalues described 
above. 
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• For 3*7 > 0, one has 3?Ao(0 > 0, ±3?(A±(C)) > 0, and for > one 
has 3?Ao(l) > 0. The matrix — Mo(^) has three eigenvalues of positive real 
part, and two eigenvalues of negative real part. 



The associated eigenvectors are given by 



Ea,+ {0=-Pi+{P + lOi2+f3i3 



(30) 



// we introduce 



( 7f \ 

P 



V y 
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\ 



-e 






the eigenvectors F± associated with A±(^) are given by 



F+ 



ii + Q^" + 



a 



/ 



■i?o(0- 



(31) 



A± - "7^ (A± - ^70^ 
The proof of this Proposition is straightforward, except for the sign of the real 
part of A±(^). For this, we use 



2v 



where ^ > (if A < one uses [A + iBY = {-A - iB^ . and A = Q leading to 
{iB)"^ = —B"^ is not possible when JR7 > 0). Hence one obtains 



A^-B-" = ^+a2^2^a^sR7, ^2 + ^2 = +a^p^ + a£,lR^f + a^e{'^lff^ 



^2v 
1 



As A±(e) = -i-±{A + iB), 5J(A+(e)A_(0) = (-^ + A){-^ - A) + B\ 
and 3?A+(0KA_(0 = 3?(A+(C)A_(C)) - = -ae^+i3?7 - a^/^a _ 52 ^ q, one 
obtains that the product 3?A+(^)3?A_ (^) is strictly negative, hence the real parts 
are of opposite sign, hence 3?A+ (^) = + A. 

This calculus also defines uniquely in the case 5i7 > the eigenvalues A+(^) 
and A_(C). 

Note also Ao(a,/3,7) and -60(0;, /?, 7) the quantities such that Ao(a,/3, 7) > 
and (Ao(a, 7) + iBo{a, jS, 7))^ = i + + aj. 
Using ^ = ^''"'"^(1 — C)) we also introduce a new system of unknowns : 



dy 



Y = TX 



^2 - ^iXi 

X3 — aPx4 
j^a(3x4 
X5 



(32) 
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and the system on Y, equivalent to ([^H) is 

^ + a(7 - |)yi - a/3y3 + a{l - Oi^ - ^ - f )^4 + a^x, = 
< ^ + f yi + «/3y2 - a7ey3 + a(l - 0(^212^)^4 + a/3x5 - 
^ + ajZi^^{y,-^x,)^Q 

namely 

dY 

— + aB{tf3,j)Y = (33) 

where 

Bii = --fii + (7 - |)i2 + fis + i-Z{»4 

< = - /3i2 - 7^13 + /3z5 

Sz4 = (1 - 0[(/3 - + (^(7 - ^) - f )*2 + (i + 7^)^31 + 7*4 

The next section is devoted to the statement of the methods used to find the 
solutions at infinity for systems of ODE which coefficients depend on y, and of 
the general set-up to find solutions bounded at ±cxd. 



2 Evans functions and application to the Kull- 
Anisimov system 

2.1 General framework 

This section recalls results of the paper of Alexander, Gardner et Jones [2], as 
well as the methods developed by K. Zumbrun D. Serre [24], S. Benzoni- 
Gavage [2j and other authors. Its purpose is to study solutions with a prescribed 
behavior at infinity of an ordinary linear system of differential equations. It is 
used in particular to identify solutions going to as y ^ ±00. 



In the general case, we consider the ordinary differential system 

^=A(t,a)2/, (34) 
when A is a regular matrix (for example analytic in a). 

We notice that the vectorial product yi A 2/2 of two solutions yi and 2/2 of (|34p 
is solution of a new differential system on h?{K'^'') which matrix is denoted by 
A^'^\ because 

-rXyi A y2) = {Ayi A + yi A Ay2). 
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Similarily A 2/2 A ... A is solution of an ordinary differential system on A*^ (R" ) 
whose matrix is denoted by A'^'^^ . The matrix A'^'^^ is given by 

fc 

A'^''\et, A e,,... A e,J = ^e.^ A .. A Ae^, A ..e,,. (35) 

1=1 

When the matrix A is diagonalizable, with eigenvalues Ai < X2 < ... < A^ then 
^(fe) 

is diagonalizable and its eigenvalues are 

iG/,/c{i,....d},Card(/)=fe 

The largest eigenvalue of A'-''^ is J2p=i Ad+i-p, and its smallest eigenvalue is 
Ap- 

Under the hypothesis that Xd-k < A^+i-fc, the largest eigenvalue of A^'^^ is 
simple. Its associated eigenvector is the vectorial product of the eigenvectors 
associated with {ed+i-k, •■, Cd). 

Recall that the space A^'^^'R'^) is of dimension 1, hence the matrix associated 
with A(^) is a number, which is equal to Tr(^(i, a)). The associated differential 
equation is 

^(2/1 A ... A yd) = Tr(A(t, a)){yi A ... A yd). 
Hence the vectorial product of d solutions of ([51)1 satisfies 

j/i A j/2 A .. A a) = j/i A 2/2 A .. A ?/d(io,a)e2:p( / Tr{A{s, a))ds). (36) 

Jto 

whose solution is the Wronskian of d solutions of the system. 



2.2 Notations for the Kull-Anisimov system 

Remark that, v^rhen 7 is complex, one has to replace R by K — C but 
nothing will change as what is important is that we study objects on 
a field, which can be IR when 7 is real, and which is C when 7 is 
complex. 

In the set-up of this paper, the matrix —Mq admits three eigenvalues of 

positive real part, which may be associated with the solutions going to when 

y goes to +00, and has two eigenvalues of negative real part, which help to 

understand the solutions going to when y goes to +00. 

It is in general hard to compute the solutions associated with an eigenvalue 

(2) 

of the matrix —Mq. However, we may compute the solution for the matrix Mq 
associated with the smallest eigenvalue Xa.- + A_(^), and the solution for the 
matrix Mq^-* associated with the largest eigenvalue Ao(^) + A-|_(^) + Xa.+- 
We introduce from now on the base vectors in A'^{K^) which take into account 
the role of in B: 
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(38) 

I 





51 = n A «2,52 = ii A 43,^3 = zi A 15 (37) 

54 = 12 A 43,55 = 12 A 15,56 = A 15 

To these vectors are associated the fohowing vectors in A'^{K^) such that fi A 
= «i A «2 A 13 A «4 A 15 = 5j A 5^, V«, j. We have 

/f'- = 12 A 13 A 15, /2" = -il A 13 A 15, fs = n A 12 A 15, = ii A 12 A 13 

5^ = ia A 14 A 15, 5^ = -12 A 14 A 15, 5^ = -12 A 13 A 14 
5^ = »i A 14 A 15, = ii A 13 A 14, 56" = -«i A 12 A 14. 

f 2 ) 

Note that we shall use in the sequel the eigenvector of the matrix A^q "^(1) 
associated with the eigenvalue —a/3 + A_(l): 

W+ = (1) A Ea^_ (1) = /3/i + (7 - P)f2 + Ph -I3.U + G (39) 

where G belongs to the space generated by gj,j = 1..6, and W+^i = (3, VF+,2 = 
7-/3, W^+,3 = /3, iy+,4 = -/3. We also introduce n{a) = -/3 + ^^i-^ such that 

A_(l)-a/3 = -l + aAf(a). (40) 

We notice that ri{a) = —B — 1 ^T" ■„ , hence 

Vae [0,«o],/3e [/3o,/3o~'],l7l </3o"\5R7>0,|/i(a)| <2ao/3o"2_^3/3o-i. (41) 
We construct a solution of the system 



+ A-4^^X(2) = (42) 



which belongs to the family of its most decreasing solutions when y goes to +00. 
Siniilarily, we construct a solution of the system 

^ + Aff X(3)^0 (43) 
dy 

belonging to the family of its most decreasing solution when y goes to —00. 
It is useful to introduce the following transformation for the study of the solution 
when y — > +00: 

y(2) ^ 2,(2)^(2) 

where 

Lemma 3 Let ^ d ^ tjij he two solutions of \2b]) . Ve denote by > = 

Vifi+V2,f2 + v^J^ + V4^Jn + Yl^j=i ''^39] ^ solution of i33]} . The associated solution 

y(2) ^ 2^(2)^(2) ^^j^g^ 



6 

^ Zi/i + Z2/2 + Z3/3 + ^4/4 + 
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with 



Ml 




a^iv2 


f Wi 


M2 


= a^vi — 




-t- W2 


M3 


= a'y£,V4 - 






M4 


= aPv2 — 




W4 


Ms 


= $^V4 + 


W5 




. Me 


= aPv4 + 


We- 





For the construction of the solutions when y — s- 
transformation of the unknowns 

^zi = xi, Z4 — aPx4. 

The system deduced from (f26| is thus 



-00, we use the fohowing 
(45) 



+ 2q!/3zi + aPx2 — 017^X3 + £,'^Z4 = 



Introduce 



^ + af3zi + C"Z4 - a(3x5 = 
^ + a/Jxa - a/32:4 = 

t — —a(3y. 



(46) 



We have 

Lemma 4 Let Z^{t,a), Z'^{t,a), Z^{t,a) be three solutions of (JS^ - Write 

4 6 

= zi(t,a)AZ2(t,a)AZ-''(t,a) =;^/,(t,a)//+^gp(t,a)g^. (47) 

T/ie solution of J^^l i associated with Z'^^^ is 

w^^\y,a) - /i(-a/3y,a)/i^ +C(y)E-=2/i(-"/32/'«)// 

The proof of these two lemmas is straightforward. 

2.3 Construction of the Evans function for the Kull-Anisimov 

system 

Recall that the Evans function is characterized by the vectorial product of five 
solutions of the system. As it writes X'-^^ A w^^\ we notice that, from Lemma 
[3] and Lemma |4l we have 

a/3ex(2) A = ^{Migi + M252 + M353) + ^^(^454 + M555 + Mege) 
+ (1 - OiZiifi - j,9i - 92) + e^2(/2 + ^51 - 9a) 
HZsifs + ^92 - ^54) + ^^4(/4 - ^93 - ^95 - ffe)]- 

(48) 
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Let Co be the limit of (1 — ^(y))e^ when y — + +00. To ensure uniqueness for 
the systems p2)) and p3)) and to adapt the constant in the system on Z,M, we 
consider the solution w^^^ of the system such that 

converges to W+§^ when y — > +00. Similarily, we consider the solution it;^^ of 



the system p3)) such that 



(50) 

Let B^^-' be given by 

' A = (7 - i)/2 + f /3 +(1 - Olii - f - ^)9l 

+(^+7^^)54] 

5(2^/3 = - /3/2 + 7(1 - 0/3 - Pk +(1 - + ^ )ff2 

+ (f + - ^ )ff4] 

I/2-/3/3 + 7/4 +(l-0[(/3-^)ff3 
+ (-f-^ + ^)55 
+ (|+ 7^)56] 
-751 + P92 - §54 
-/351 -7(l + 0ff2+/353 

+(7-^)54 

- + fi) +f.9i +/552 -753 

+(7- ^)55 + fg6 

- 7C54 + /355 
-7C51 + P9i + /356 
-7C52 + /3C53 - 1^4 
-/355 - 7^56- 

(51) 

c,S(2)y(2) ^ (52) 



S(^V4 = 


-7^1 - 




l-fJ2 




13 f 
-135/3 


5(2^53 = 








B^^^g, = 


-1^/2 




-1^/3 



The system on y^^^ is 



The Evans function that we shall use is given by: 
Definition 2 Introduce 

Ev{a,l3,-i) = a/3w^^'(y) A ii;L^'(y)e-/o'("'^«2'')-(«(2''»")'*y'. 
This function, independant of y, is called the Evans function of the problem. 
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It is easy to derive the 

Proposition 2 The complex number 7 is an instability growth rate according 
to Definition [7] if and only if 

Ev{a,l3,-f)^0. 



The proof of Proposition [5] is to be found in [5] . 

Reduction of the Evans function We deduce from the systems and 
that w = w^^^ A is solution of From TrAfo = ~aj^ + T = 

-ajS, + I + wc deduce that the derivative of w(2/)e- /o'("'rC(y')-(«(a'))'')rfy' 
is zero, hence 



Using we obtain 



w(2/)e- Jo "-r^t!/ constant. 



£^w(a,/3,7) [Ziifi - j^gi -52) + ^^2(/2 + ^gi - 34) 

+?^3(/3 + ^92 - j^gi) + iZi{fi - ^gs - -^55 - 56) 
I <r Migi+M2g2+A^3g3 

I c2 Migi+M5g5+Mege ]„- f" ay£(v')dv' l-C(O) 

1-4 J'^ e(o) • 

2.4 Statement of the principal tools for the study of the 
Evans function 

The aim of this paper is to compute the Evans function through the calculatiord 
of and . 

Theorem 1 • For all > 0, there exists ao > 0, /3o > such that, for a < 
ao,P € [/3o,^],7 € ^gj;] ^"^^ U such that ^{y) G [Co,l] there exists a unique 

solution li;!^^^ of satisfying ^4^ - 
• Let 

w = r(2)?«f (y) exp(-a^(a)y) - + W^+,2/2 + + VF+,4/4). (54) 

T/ie function w is analytic in {y,a) for ^(y) £ [^q, 1[ and for a < uq. For all 
ao > and for all ^0 g]0, 1[ there exists a constant C{^o,ao) such that 

yy,m G]Co, i[, \My)\ < c(Co, ao)(i - ay))- 



^Note that the techniques of diflferential equations allow us only to compute the solution 
of a differential system with non constant coefficients only when this solution is associated 
with the largest eigenvalue of the matrix in the neighborhood of —00 and with the smallest 
eigenvalue of the matrix in the neighborhood of +00. 
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• There exists ai < ao and R > 0, depending only on ^q, such that for C^o < , 
w admits an analytic extension for y such that y G [(^)^jCo]- 

• For ao > and (3o > given, and for all a < ao, (3 S [(3o, ;gj^],7 G [0, 

there exists tg > such that there exists a unique solution w;^^^ (y, a) of for 
—aPy G [<o,+oo[ satisfying i50\) . Moreover we have 

w^'\-^,a)exp{~ f ^\\++a{(3 + -fay')))dy')-Wo{t)e^H-^ ^0{ai) 



afi 

uniformly for t G [ip, 



(55) 



We see in this theorem the division in three regions for the computation of the 
solutions of (|42| and of (l43l) . The first region ^ (y) g [^o , 1 [ is the aim of Section 
[21 the second region (which extends the result of [^qj 1[) is studied in Section [231 
and the solution in the neighborhood of — cx) is characterized in Section [5] In 
the next paragraph, we describe the systems that we shall use in what follows. 
One of the main problems is that the problem to solve is not a Cauchy problem, 
but we have informations at ±oo for the solution of (j26[) that we want to study. 
Let us say a word on this system. As ^ goes to when y goes to -co, we notice 
that the matrix Mq is singular when y goes to — cxd. 

Moreover, the term ^ prevents us to have a result which is uniform when (3 goes 

to 0. A possible choice to overcome this difficulty is the choice /3 G [/3o, 
Finally, even if we remove the singularity of Mq at j/ — s- — cx) by whatever 
method, another problem is induced by the behavior of when y goes to —oo, 
because ^ goes to as This will lead to a fuchsian problem (Section [5]) . 

Note that in the theorem [T] we defined unique solutions of the problems (|42l) 
and ()43p with the prescribed behavior at infinity. 

The behavior of the solutions induced by the theorem [T] lead to another expres- 
sion of the Evans function of relation ([55]) . 

Assume that y < and denote hy t = —a(3y. We introduce the functions 
Lp{t, a), Rj{t, a) given by the equalities 

5p(t, a)e'*+^ ^4 = Lp(t, a), 1 < p < 6 ^^^^ 

/j(t,a)e^*+*^*o«(-^)'^*'r577 =i?j(t,a),l < j < 4. 

The corrective factor e^*^^ -^^o ^^"'^^'^^ t^^ is induced by the relation (|55|) . 
Similarily, we introduce 

z,iy,a) = Z,(y,a)e-"^(")^mp(y,a) = Afp(y, a)e-"^(")^ (57) 
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We obtain the relation 



to 



Ev{a, (3, 7) = e-(2+^)*t^ exp(- ^ a^S,{y')dy') 

- ^Li - L2) + ^Z2{R2 + ^Ll - i4) 

+fz3(-R3 + ^^2 - ^^4) + e^4(-R4 - ^^3 - f ^5 " ^e) (58) 
I ^ mi_Li+in2-L2+'»3-L3 

I <;2 m4L4+m5L5+m6-t-6 1 

1"C J C(0) • 

In what foUows we shall describe the functions R, L, and z, m. 



3 Calculus of the solution in the overdense re- 
gion 

The aim of this section is to obtain w^_^_ ' , which is the unique solution of the 
system ([^^ under the condition (HS)) in the region [?/07+oo[ for all t/o- The 

first idea would be to try to apply the Gap Lemma. However, the eigenvalue of 

(2) 

smallest real partof — Mq (1) is Xmin — ^^-(1) ^ ctP^ sjid the next eigenvalue 
is A* = A_(l) + aj. As A* — Xmin — Oi{l + P), the difference between two 

(2) 

eigenvalues of Mq "^(1) is not uniformly bounded below for a e]0,ao[. Hence 
the hypotheseses of the Gap lemma theorem are not fulfilled. 



Remark As the coefficients of the differential system behave as 1 — e ^ when 
y goes to +00, one may obtain K^^^ through a Volterra expansion as 

r(2) (y) = yle(^- e-yr{y, a, (3, 7) 

r being a remainder term. We prove that the result that we obtain is analytic 

(2) 

in a for /3, 7 in a certain compact set. Once we know that, we identify w\_ 
(deduced from Y^^^ and obtained also through a Volterra expansion) and we 
express it with an expansion in powers of a, which is bounded by a geometric 
series for £_{y) e [^o,l]- 

We then plug this developement in the differential system and we identify the 
terms. We check that the coefficient of in the expansion in a of the solution is 
of the form ^„}+'i Xj{£^), the radius of convergence of the series depends on ^(yo)- 
The aim of the explicit calculus is to deduce a new region of convergence of the 
series using the behavior in g"+2 ^j(C)- In this new region of convergence, we 
have a converging series, which coincides with the original one for ^ G [S.{yo, ![• 
Hence it is the extension of our solution in the new region of convergence. This 
is an important feature, because it helps to have an overlapping region between 
the most decreasing solution at —00 and the most decreasing solution at +00. 
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We introduce Zj and mp through 

j=i p=i 

where /i(a) has been given in (|40|) . We recaU that Zj W+j and rrij when 
y +00. 

The following Theorem summarizes the results of this section. Let us introduce 
eo >0. 

Theorem 2 1. The functions Zj and rrij have a normally convergent expan- 
sion in powers of a of the form: 

Zk{y, a) = E,>i e^flfej (0 + W+^k,k = 2, 3, 4 
■mi(y,a) = Y.j>i j^bij{£_),l = 1,2,3 
mp{y,a) = Y.j>i ^^bp,AO,P = 4,5,6. 

We introduce Si — 0, 62 ^ 63 ~ 64 = di = d2 = d^ = 1, d^ = d^ = de = 2. 
Let K be a compact subset of M*^ x R. There exists R > and ai(^o) 
such that for /3, 7 in K and a < ai (^0) we have for ^ € [^Qj 1] 

I d l bi.,{e,) \\ ^ w 

\hiAO\<R'{^-0- 

Assume this theorem is proven. We show that, for all there exists ao(Co) 
such that, for a < ao('Co) the power series is convergent for ^ G [^0, !]■ Consider 
now < 1 and < 1. The power series J2j>i ^-^fUp defines an 

analytic function which is the analytic extension of the sum of the normally 
convergent series J2j>i ^<^7+dc for < ai(^o) and ^ G [Co,l]- A similar 

result holds for the series J2j>i ^Jj+ii ■ 

The proof of item 1) of the theorem is in Annex [71 We calculate in the first 
subsection of the present section the first order terms of z and m. 
In the second subsection, we prove by recurrence the structure of the j— th term 
of the expansion of z and m, which is a consequence of the structure of the 
system. We deduce the precise estimates wich help us to extend the expansion. 

3.1 A new formulation of the system 

We introduce Zj,l < j < 4 and nip , 1 < p < 6 given by 
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a/3 



4 6 4 6 

-T^^)wf{y,a) = (y,a) = ^ Af,5p = e^^^-'^E ^^-Z^+E "^p^p)- 

(59) 

The relation w^^^e'^e-^''^")^ ^ F_(l) A Ea-{1) = W+ imply that Zj ^ W+^j 
because (■]^_|-)gy — *■ 1 when y —f +oo. The system on Z, M writes 



f + fJoZ + aBiOZ + ^(Lo + C^i)M = 
dM ^ ^j^^M + aD{i)M + ai^C(OZ = 0. 



(60) 



that is 



/3CZ3-7e^4-^ = 

(7 - ^)^i+ 7^2-/3^3- 1^4 

f + /3Z, + 7(1 - 0^3 - /^^4 ~ P^'YJ"'^ ^ 

-/3Z3 + 7^4 + f4|=0 

i ^ -7M1 - /3M2 + |M3 + ^ (-/3M4 - 7M5) 

+ (l-0[(^-]^-f)^i + ^^2]=0 
i ^ +pM^ _ ^(1 + ^)Af2 + /3M3 - 7P/r6 

+ (l-e)[(^+7V)^i + '^^3]=0 
li^M, +^^^^_^M3 + /3CM6 + (l-a^^4 = 

IdM, _|^^^ _^ _ 1 _ ^^^^^ ^ ^^^^ _ 1^^^ 

+ (1 - e)[Z2(i + 7ij^) + ^3(f + ^ - ^ )] = 
+(7 - ^)M3 + /3M4 - /3M6 + (1 - 0^4(-f + i(7 - ^)) = 



idM, ^1^^ ^ _ ^^^^^ ^ _ ^)^^( 1 ^ ^ 0^ 



(61) 

This system writes 

-'^+PK{r,^{y),P)U ^0 (62) 
a ay 

where — (Zi, Z2, Z3, Z4, mi, TO2, m3, TO4, ms, mg). The system on (z,m) is 
deduced from ([^ by replacing the matrix (3K{r,^{y)) by (3K{r,£^{y)) + /x(a)/. 
We verify that — > when y — > +c». In what follows, we describe the analytic 
expansion of the solution in a when a is in a neighborhood of 0. 

3.2 First terms of the expansion in a 

From 

Zj ^W+j,mj ^ (63) 

we obtain z° = W+j , mJJ = 0, that is zJ* = /3, 2° = 7 - -^3 = /3, 4 = We 
replace in the system these relations to obtain the system on the term zj and 
mj. We get 
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If + m(o)z§ + (7 - j^yi + l4 P4 H = 

M + ^(0)4 + + (}-4 + 7(1 - - fizl = 
^+M(0)^4°-/?^3°+ 7^2 = 

^ + (i-e)[(^+7ii^>? + ^^§] = o 

+ (1 - + 7 V) + ^3(f + ^ - ^ )] = 

i^ + (l-0z2(-f + ^(7-^))=0 
As /z(0) = — /3 — 7 we deduce 

# + ?-i = o 

1^ + ^ + 7/3(1-0=0 
^ = 

^ # + (l-O[(^-^-f)/3+^(7-/3)]=0 

# + (l-O[(^+7V)/5+^/3]=0 
^-(1-0^/3 = 

^ + (l-0[(7-/3)(^+7ii^)+/3(f + ^-^)] = 
i5i-/3(l-O(-f + i(7-^)) = 
_ ^-/3(l-0(i+7ii^) = 

We obtain 

r zi = /?(/? + 7)^ 

Ui = o 

and the following system on mj : 

^ + e^K + 7/3(l-O + (7'-/3')C]=0 

^ ^-e^(/3^-7^) = o 

^ + e^[(7 - /3)(| + 7(1 - 0) + (/3' + - 7^0] = 

^-e^(-/3^+7e(7-^)) = o 
. ^-e^(€ + 7/3(i-0) = o 
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Hence we obtain the expansion of the solution for ^ > 
r ^i=/3-^(l-r)/3(/3 + 7) + 0(a2) 

24 = _/3 + 0(a2) 

The next subsection is dedicated to the precise study of the behavior of w\_ , 
which depends on inverse powers of ^ and cannot be extended directly to ^ ^ 0. 



3.3 Uniform estimates of the solution w^^ for ^ G 1] 

We write 



So. 



(64) 



We introduce the constant matrices Jo, and such that J — Jo+^B{£^), L 
Lo + ^Li,K = Ko + £.D{^). The system ® yields 



^ + an{a)m + fKm 







+ a^^Cz = 



(65) 



The aim of this paragraph is to find a simpler formulation for the unique 
solution of this system. 

As the solution going to at infinity of ^ = is f{y) = 7;qri(l — 
direct estimates of the behavior in ^ of a coefficient Uj or rrij of (|65p lead to a 
multiplying factor of the form 

a 



v+l ■ 



The behavior of a^'^^Uj+i or of a^+^TOj+i in ^ is then given by for the 

next coefficient. However, the structure of the system allows us to obtain a 
lower inverse power of in the expansion in a. For this purpose, we introduce 
the new unknowns Op j , bq^j such that 

z^ = W+,^ + f^^^AO + z^+\a,0 

zp = W+^p + Ef^i ^ap,,(0 + z^+Ha, 0,P = 2, 3, 4 
= Ef=i J^hAO + m^+\a, e), fc = 1, 2, 3 



(66) 



N 



,N+1, 



a,0,^ = 4,5,6 



The previous quick analysis would suggest that Opj+i is of order i when every 
ak,j,bij is bounded when ^ goes to 0. This is not the case, and the crucial 
equality states as follows. We introduce the diagonal matrix T such that 
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T 



/ 


ai 






/ 




\ 




0-2 














03 














a4 














bi 














b2 














b3 














bi 


























v 


&6 


) 




V 




/ 



There exists a matrix C(^,/3,7), polynomial in ^, such that 

T = TC. 




(67) 



We shall make use of the following fundamental Lemma, noting that at each 
step we solve an equation of the form 



dy 



where a = vj + d, d ~ 0, 1, 2. 

Lemma 5 The unique solution going to when ^ goes to 1 of 

df _ 



dy 



^■s I{y) — drj. We have the 



estimate 



|/(y)|<^||A|| 



Proof From 



\fiy)\ < 



.^+1+0 



—(—(1 - 

^ d ^ 1 

J d^ ly + a 



p^)), we deduce 
1 



(1- 



iy-\-a 



))rf^=ll^l|c 



1-C 



The equality Sl^Z-i = + q,)(i + _ \)Y+'^-^dt implies 



1| < {y + a){l-0- 



(68) 



hence the lemma. The indices defined in Theorem [5] will express the weight of 
each coordinate of the vector U defined in 



1 
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In the system ([M)) . write 



ml \ I ^ P^i \ h. 



we get, for the term in a^'^^: 



rfyV^<.0 + l) \^ hj + l J' ' 1 \ Q J Z-.^ y bl J ^-^J \ hj 

This system of equations becomes, for j > 2: 

_±( _apj±2_\ I iI2l=ot^i"-P.j-iC^) + {Ciiaj+Ci2j^))p _ 
l^ap + .-(j + i) J "I ^t,p + uj r Hj-lVy + ^p — U 

d I bpj + i N (II]i=o^i^p,j-i)+C2l(l-^)aj+C22bj) _ 

and the equahty for j ~ I 

d( Op, 3 + 1 X (Z;i=oA'!ap,i-i?"') + (Ciiai+Ci2T-^))p ^ .^ , _ 

d / bp,j+l ^ I (Ei = o^i''p,j-')+C2l(l-C)aj+g22fcj) , ^1 _n 
d^l^dp + .-0 + i) ; ^ ^dp+i-O+i) r 1,1 ?j'^4+p — u. 

We have the identity Ci2(l)&^(l) = + fij^iW+j, which is necessary to 

obtain that the source term in the equation on flpj+i vanishes at ^ = 1. 

3.4 Behavior of the terms of the expansion 

The regularity of the quantities ctp.j, bqj is given by the foUowing proposition, 
which gives precise estimates on the functions provided that /3 and 7 stay in a 
compact set: 

Proposition 3 Assume that /3, 7 are in a compact set K of R+ x C , namely 

Po<P<^Al\ < -^,3*7 >0. 
Po Po 

Assume that S^q g]0, 1[ is given and that there exists ao(^o) such that the analytic 
expansion of the solution of \65\) is valid for a < ao(^o) ond S^{y) e [Coi 1[- 

There exists R and M depending only of ao and f3o such that, forall p = 
1, 2, 3, 4, forall q = 1, 6, for all j > I we have the estimates 

If we introduce aj — aj{0) and f3j = bj{0) we have 

K(0-a,(l"OI <M^C(l-0 
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From this proposition we deduce the 

Proposition 4 Under the same hypothesis as the proposition\^ 
i) there exists q;i(^o) such that, for a < ai(^o); the functions 



are analytic for < . 
2) For y such that 



the functions 



Up{y,a) ^W+^p + ^ ^''Ap{j^j^,a),Vq{y,a) ^ ( ''^ Bq{j^^--^ , a) 



are solution of the system i65\) and extend the gap lemma solution. 
3) Introduce 



The functions Ap and Bq are analytic for C, < , and we have the inequalities 
forC<jj 



\BqiC,a)\<{ 



C 



C 1 - MC 



l-MC' 



We deduce, for a < ai{^o) and C < ig: the equahties 

I aiUpiC, a) = C; ET=i ^P,i(' + «; ^P(C, a),P = 2, 3, 4 
a^VqiCa) = C- J2JLi f^uC^ + oi"Sq{C,,a),q = 1,2,3 
^ a^VqiCa) = C~ T.'jLi l^idC' + oi"Sq{C,,a),q = 4,5,6 
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Proof of Proposition [4] Let €]0, 1[ be given. There exists ao(Co) such 
that, for a < ao{£,o), the solution satisfies the conclusions of the gap lemma 
(which means that the solution of the system (|65p is analytic in the region 
a<aoi^o) foree Ko,!])- 

Moreover, from proposition [31 for all there exists R depending only on /3o 
such that for C < the functions Ap{(,a) and the functions Bq{(,a) given in 
Proposition |4] are analytic through their expansion in a for C < 

Introduce ai(Co) = min(ao(Co), For a < ai(^o) and C < ^, ^ = (|)- > Co- 
This means that 

(^+.P + (-)~^p((7)",«),(-)~S,((-)^a)) 

is solution of the system ((65)) for a < ao(Co) when y is given by ^(j/) — by 
construction of the analytic solution given by the gap lemma. 

By uniqueness of the solution which is analytic in a, we check that this func- 
tion is also solution of the system ([65]) for a < ai(Co) and ^{y) e [{aR)~ ,S,o], 
because the analytic expansion defining the solution in the set-up of the gap 
lemma can be rearranged and the remainder term is regular enough (and uni- 
formly bounded), and because the two solutions are equal at the point y such 

that ay) = 

Hence we extended the solution for in the interval [{aR)~ , 1]. Proposition 
3] is thus a consequence of Proposition [31 



Proof of Proposition [3l We prove by recurrence the inequalities (TfT]) . (|72l) . 
((73|) below: 

I " ( "i-:> )\ < tl (70) 

We assume that (/3,7) K d [/3o, ^] x {7, 17I < ;9(7\3?7 > 0}. 
First step: From the identity 



and from the same identity on bj, relying on the fact that aj and bj are at 
C = 1 we deduce the inequality 
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We use to obtain 

l«pj(0 



or 



I MO, 



< i?^ (74) 

< (75) 



and of course 



1^(1)1 -||(|fe)(l)l<^^- (76) 

We assume that ([7T|) . ([72)) . ((75|) are true for / < j. We deduce the inequahties 
dZlD, dZSl), (El]). 

The equation on Opj+i rewrites 

+(Ci2^ + Ci2(l)6;(l))]. 
Recall that we have the equality 

7 + a/32 



1 + \/l + 4a(7 + a/32) 



hence for (/3,7) in the compact the function /i(a) admits a DSE at a = 0, 
of radius of convergence greater than = min(l, ^). Moreover, denoting by 
Co = Y.Zo lMi|6'o7 we have 

CSO 

1=0 

Using the norm of the matrices Cn + fiol, C12, C21 and C22 + Mo^ (with < 
^ < 1) and the inequalities ([7^ . ([75)) . ([7^1) we obtain 



< ^wTTTTTTwELi A^^i?-'-'^' + (3|7l + 6/3+ + 
< Twrrfc^fCo + 3|7| + 6/3 + ^ + ^] 



as soon as 

<min(l,^). (77) 



We have the same type of estimates for bqj- 



^ ^.o+D+i+d, iCo + 3|7l + 6/3+ ^ + 
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Hence there exists a constant D(/3o) depending only on the compact set K, such 
that for R satisfying (|77p and assuming (ITTI) . ((72|) . ((73|) at the order j, we obtain 
the estimates 

I J.'V + n ''I — 



l^f Vd„+l/f7"+l) ''I - 



The last estimate that we need is based on the difference of derivatives and we 
use the identity 



^^^^ + /'(I) = ^ /'(/'(I) - /'(I + s(e - mds 4^(1 - D 



for /(I) = and / e C^. Assume the inequalities - /'(1)| < C'fe^ and 



|/'(1)| < C. We obtain 



We apply this estimate for /(^) = ^Sp+^i and a = u + 1 {(3 = dp + v{j — /)) to 
obtain the inequality 



We thus deduce the estimate 



I d , bqj+i _ d bqj+i ^ RW{/3o){l-^) 

We use the estimate for j = 1, for which there exists a constant Di{Po) such 
that 

I d flpj ^ -Pi(/3o) I ^ -Pi(/3o) 

._&g^w.s _ d bq^i ■Pl(/3o)(l-0 

It is then enough to consider 

R^miixil,^,D{po),Di{l3o)) 
Po 

to obtain the inequalities ([7T|) . ([7^ . ([75]) for all j. We write each term Upj = 
'^p,j + ^Cpj, = l3q,j + £,Cqj, and we have similar inequalities for the terms 
Cpj and dqj, with a coefficient M depending only on the compact set K. To 
obtain the estimate of the rest (versus the leading order term), we denote by 
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Cj and dj the functions such that aj^p{^) = (1 — C)(Q^j,p + ^qjiO — 

(1 — OiPj.q + ^dq_j{C))- We prove in a similar fashion the inequahties 

Irf^^^s;:?^^! - ^^OT ^^^^ 
l|(^)«) + 4„(i)l<^^ ,S0) 

This ends the proof of Proposition [3l 

4 The simplest discontinuity model 

Before studying the couphng between the hypergeometric region and the over- 
dense region, we shah in this section study a simple model where the profile of 
density is 

(,{yo),y e] - oo,yo] 
C(y),y e]yo,+oo[- 



It is a slightly better model than the discontinuity model (see [21j) for two 
reasons: 

i) we assume that the density profile is continuous, 

ii) it corresponds to a simple form of the energy equation. 
The stationary associated quantities p* and u* are solution of: 

= -p*{x)9 

-£{uoix) - L^Va±{Z{p,{x)))) = -K^'(PaCo)Par(2/o)'5.-,oLo 

coming from 

dtp + d\Y{pu) = 

dt (pu) + div(pM ® u + pid) = pg (81) 

div(u + LoKVZ(p)) = Va^S.,^y,Lo- 

It is easy to see that the stationary solution is thus given by 

1 -c|),ye]yo,+oo[ 



Poiy) 

We have 



- Pag^oLoiy- yo),y e] - oo,yo] 
'j^ ~ Pa9 Jy^as)ds,y e [yo,+oo[ 
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Proposition 5 The linear growth rate associated with the system \81\) in the 
case Ho fixed independant of a in the regime a — > and (3, 7 fixed is 



1 = 

which corresponds to 



l+Co Co 



f^P^_PO _ kVblowoff- 
Pa+PO 



We recall that for all S^c- 



Note that the leading order term in a conies from the coefficient of 14 and of 
T". As the leading order term in a of the solution in the region [j/o,+oo[ in 
the case where j/o is independant of a is given by the leading order term of 
/_(1) A i?a,_(l), we find 14 A £'a,_(l) = W+ as leading order term. The Evans 
function is 



Ev{Q, Co) = W - 7Co) .2 ^°2.2 [/3(Co - l)iT^ + T°) + {W + 7(Co + 1))T{'] 

P 7 ?o 



hence 

i?«(0,/?,7,Co) = /3(/3-7Co) %2 ^° 2.2 [Cg(l - Co) - (1 + Co)(/? + 7Co)'^ 

P - 7 Co 

The theorem [5] is proven. 



Remark This result is somewhat surprising: in fact when Co ^ the limit of 
the growth rate disappears because it becomes negative: in the previous equality 
Co must verify 

Co' + - l)Co + > 

to obtain a positive growth rate. It is easy to check that this is never satisfied 
if > max^gg]o^i](CoT^|^) and is satisfied only in a region C- (/?), C+C/?) for 

/32 <max^,g]04](Co'i+|)- 

We may thus conclude that this model is not relevant according to the phys- 
ical results. 
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5 Solution in the hyper geometric region 



In this section, we identify the solution when y goes to — oo. We show a result for 
t < to finite. As the coefficients behave as |y|~ " , we call this problem a Fuchsian 
problem, and, as the leading term of the equation leads to a hypergeometrical 
equation, we call this the hypergeometrical set-up. We introduce 



.Pit) 



(82) 



Let 



k=l 



k 7 



(83) 



Note that h is a complex valued function because p is complex, but as ^ and 
jh{p) are real, ph{p) is also real. From lemma [Jl we deduce 



[i+pKp)). 



(84) 



The system of equations on Z 



( zi\ 

\X5 J 



obtained from (l46l) is 



dt 

dX2 

dt 

dX3 

dt 



\{l - al-q)zi 



P 
-pzi 

- 2zi 



X3 



"+^7.4, = 



X2 - PX3 



+ zi + ^2:4 - = 



^Z4 = 



dt 

dx; 

dt 



X3 — Z4 = 0. 



(85) 



Introduce the matrices 



Moip) = 



/ 





1 








\ 




/ 


Zl 


\ 


p 





-1 



















2 


1 


-p 










,NZ = 




Z4 




1 











-1 








Z4 




I 





1 


-1 





) 




\ 





/ 



Using the equalities 



f = ^,+j;^pih{p){l+phip))i 

ai^-^^-^{i+phip)y+ipi-^ 
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there exists a regular matrix {p, P, 7) such that ([55]) rewrites 

+ Mo{p)Z + -NZ + -^—^M{p, p, ^)Z = 0. 
at vt {vty+- 

It is then natural to introduce 

We associate with the system ((85|) the model system 

f .Afo(0)Z+lA.Z (86) 



and the model system of the vectorial product of three solutions of ((86)) 

We need the following 

Lemma 6 The matrix Mq{p) lias three eigenvalues — 1,— 1. The eigenvector 
asociated with Ao(p) — —p is ea{p) = (1, — 2, — p, 0, 1). The eigenvectors as- 
sociated with the eigenvalue of multiplicity ^ A_ = —1 are e_(p) = (1,-1 — 
p, —1, 0, 1) and /_ = ^(0, 0, 0, 1, 1). The eigenvectors associated with the eigen- 
value of multiplicity 2 X+ = 1 are e+{p) — 1, 1, 0, 1) and /+ — i(0, 0, 0, 1, — 1). 

We notice that 

Neoip) = Ne+{p) = iVe_(p) = e+ip) + e_(p) - eoip) + f+ - f- 

and 

7V/+ = 7V/_ - i(z3 + *4) = lie-ip) - e+ip)) + !(/+ + /_). 
We also introduce the function 

Mt) = / V{s)ds. (88) 

P J to 

which is the integral of p and corresponds to the eigenvalue of largest real part 
of — . This eigenvalue is associated with the eigenvector P{p) such that 
eo{p) A e_(p) A /_ = ^^P{p) and we find 

P{p) = (2 + p)f^ + + - si - si - (2 + p)gi - gi - gi . (89) 

The aim of this section is to prove the 

Theorem 3 There exists a unique solution Uo{t,a) of the system I185\) and 
to > such that there exists C > such that, for all t > to we have 

|t/o(i,a)e2*+'^°(*'")| <Ct^. 
limt^+^Uoit, a)e2*+V'o(*-")t-2^ ^ p(o), (90) 
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the vector P(0) being given by i89\) . Note that this condition is equivalent to 
(Ml- 

2) There exists a unique solution of the system ^87\ l such that 

Z{t)e^H~^ P(0) 

when t —^ +00. 

3) For all to > 0, there exists C{tf)) > such that for t> to the estimate 

\Z{t)e^H-^ - Uoit, a)e2*+'^"(*'")t-^ | < C{to)ai . (91) 

In a first paragraph, me make a reduction of the system to simplify its 
resolution. 

5.1 Reduction of the system 

The system ((85|) rewrites 

— ^ Mo{p)Z + —NZ + ait-^-iR{p)Z (92) 

where Mq (p) are analytic functions of p for p < po and is a constant matrix. 
The system satisfied by the vectorial product of three solutions of (|85|) is 

^ = Mo(p)(3)z(3) + liv(3)^(3) + ait-i-ii?(p) (3)^(3). (93) 

to which we associate its model system (|87|) . 

Introduce the eigenvalue of smallest real part of Afo(p)^^\ (which is Xi{p) = 
—2 — p). Consider the matrix 

A(p)=Mo(p)(3)-Ai(p)/ 
as well as the unknown Y*(i) = yf'^^e \Ye obtain the system 

^=A(p)n + -^iV(3)n + airi-^i?(p)(3)y, (GS) (94) 
associated with the model system 

^(fe^*) = A(0)(fe2*) + -^^Ar(3)(^g2t) ^^^^.^ ^gg^ 

The aim of what follows is to identify the family of solutions of (GS) such that 
Y^,{t) ~ ct^^Eo when t — > +00, where M — and Eq is the eigenvector asso- 
ciated with Ai(p). 



39 



Remark The result that we obtain here depends heavily on the fact that the 
non zero eigenvalues of A(p) denoted by Xi{p) satisfy 

(Int)^^ / \i(p)ds ^ +oo,t +00. 

Jto 

The first transformation uses A{p) = {P{p))~^ D{p)P{p), where P{p) is 
a transfer matrix and D{p) is the matrix of eigenvalues of X{p), such that 
Dii{p) — Xi{p), Xi{p) < Ai+i(p) and all the eigenvalues are positive. We in- 
troduce U{t,a) = P{p)Y^. The system is 



^ = D{p)U + -P{p)[N(^^ + K7)^i?(p)^'^ - P{p)-^^{vA]P{p)-^U, 
at vt t dp dt 

and using the relation 



we end-up with the system 



^^Dip)U + ^^Mit,a)U. (96) 



where the matrix M writes 

M{t,a) = iV^^) +R{t,a)a-t-i, 
with the following estimate on R: 



3a°,ro > 0,V/3 e [/3o, -],V7,5R7 > 0, |7l < /3o^Mi?(i,«)l < ^^(/3o). 
Po 

On each eigenspace of D{p) we assume that N is diagonal. 

We denote by Ei the eigenspaces of D{p), 1 < i < m, Ei and being of 
dimension 1. The unknowns U are written Ui G i^^dim^i^ 

The aim of the next paragraph is to construct iteratively the solution of the 
system (l96l) . We use the methods of Levinson [18] and Hartmann [11] . 



5.2 Formal solution of the system 

It is necessary to begin with the computation of [/,„, solution of 

^ = A„ (p) [/„, + ^ (iV„„„ + ait-i R,nm ) f/m 
We consider the differential equation 

= Xm{p)Um + — {Nmm + a^t^~ Rmm)Um - f 
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As we want to obtain a bounded solution of (|96p . if we introduce (j>m{t) = 
Itf) ^rn{p{s))ds, this differential equation becomes 

If the function [/™(i)e-"*'"(*)i — ^ goes to a non zero finite limit when t goes to 
+00, then Um{t) goes to infinity when t goes to infinity under the sufficient 
condition ^"^"^^^-^ +00 for t +00, which is contradictory with the fact 
that we seek a bounded solution. Hence it is necessary (but not sufficient) that 
Umit)e~'^"^''*H ^ goes to when t goes to infinity. 
The system rewrites 

dt 

We introduce the operator 

(97) 

We verify 

I (t(") (/)) = A™ (p)T('") (/)) + ^ (A^™™ + a^i- * i?™™)T('") (/)) - /(O . (98) 
Hence the equation on Um leads to the necessary relation 

^ m— 1 

C/„(t) = -T('")[-^M„,C/,] 
i=i 

which rewrites 

-J m— 1 

t/„(i) = -^7;(")(C/,.) (99) 

with 

ii^'"^(;7,) = -^^tr('")(i-Af„,t/,). 

We replace this equality in the system satisfied by {Uj)i<j<m- We obtain 

l<J<m-l ^ = A,(p) + -5]MW(C/fe) (100) 

fc=i 

where 

M^jkiUk) - M,fe(C/fe) + ^M,„,r^™Vfe). (101) 
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Of course, we notice that 

M'^kiUk) - n.kUk - {Mjk - N,k)Uk + —M,^Ti"'\Uk) = Oia^T^) (102) 

hence its contribution is a regularizing operator. 

The scheme of the proof is the same for ah the terms of the vector U . At 
each stage, we obtain the system 



rlTT 1 '"^"^ 

1 < J- < m - e ^ = A, {p)U, + - ^ {U,) (103) 

where 



dt ■'^ ' vt 

k=l 



The operator ^('"-e) of the foUowing step is given by the solution r(™-'=)(/) 
going to at +oo of the equation 

^ = K.^MU..-e{t) + ^Ml:-ll_Mn^e){t) - fit) 

and the operators T^™ l<k<jn — e^l are given by 



t^'HUk) - -^tr(™-)[lAfil,,(C/,)]. (105) 



k . 

The construction of the formal solution is done. We end up with the remain- 
ing equation on Ui , to which we cannot apply the previous method because the 
associated eigenvalue is 0. This system writes 

j^{t-^u,{t)) = ^r^{M[r''^N,,)u,. (106) 

We deduce that there exists a constant A such that 



Ui{t) - At^ = -t^ f °° —s-^{m["'~^'^ - Nii)iUi){s)ds. (107) 
Jt 

5.3 Proof of the convergence of the previous Volterra se- 
ries 

The resolution ends up with the construction of the solution of (I107|) . For this 
construction, it is necessary to study the regularity of all the operators 
(and of all the induced operators T^^"* '^'^ and m'^'^). 

For a given function ip, such that il^it) is increasing, going to +00 at +00, 
we introduce 

A^(io) = C-([to,+oo[),3C, 1/(01 < Ci^e'^^*)}. 
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We say that e L^_g{to) if we have 

^ e C~([to,+oo[), V'(i)-5i</>m-eW < -eo < 0,t> to,t'^\'tjj" {t)\ bounded on [to,+oo[. 
These notations being introduced, we consider the operator 

K^'\f) = r -s-'^iMir'' - N,,)if)is)ds. 

Jt 

The equation that we intend to solve is 

U,-At^=t'^K('){U,). 

We notice that i7js"^(m}™"^^ - N^i){f) e L^{[U,+oo[) as soon as / G 
A Nil (*o)- Moreover, for given, there exists a constant Co such that 

\l-s-^{M[r^^ -N^^){f) < Coa^s-i-^max;eK+oo[(|/(0«-^)- 
Hence we get the inequality 

/ + 00 
s-^-^ds = Cova-t-^. 

Assume a < ao given. There exists a value of t^, given by 
such that for t > t° we have 

Hence, by induction, we get that 

|(ifW)W(t^)|<ii^,i>tO (108) 
from which we deduce the convergence of the scries j:{K^^))(i\t^) and its 

Wii 

bound by 2t " . More precisely, we have, for all a < ao 

\{K^'Y\t'^)\<^i-yt'^,t>t^ (109) 
2' ao 

hence the behavior when a 0. 
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6 The instability growth rate 



This section reUes on the relation ([58| that we obtained in the third section. 
The scope of the present section is to derive a limit, when a — s- of the Evans 
function Ev{a, (3,^). As the right hand side of ([55]) depends on t, and the left 
hand side of (|58p is independant of t, we will study, for a given (suitable) t > 0, 
the limit when a — > of the right hand side of (|58p . This corresponds to the 
calculus of the limit when a — > of the functions zi, S,{y)zp {p = 2, 3, 4), £,{y)mi, 
{I — 1, 2, 3), £^^{y)mi, Z = 4, 5, 6 and Rj{t, a), Lk{t, a) for y — Recall that 

we restrict ourselves to the regime a — > 0, /3, 7 in a fixed compact set. 
In all what follows, we introduce 

r=l (110) 

P 

6.1 Calculus in the overlapping region 

Introduce > given such that 

vt^, 1 

We introduce C,{t,a) = C(C(-^),a) = (^(^(^j^))^ - We have 

The equation on ^ gives ( + 0{ai) = hence 

C(i,0) = ^. (Ill) 

Hence we get (and it is the same for all other quantities) 

t 

a/3 

We have, moreover 



lim„^o;2i(-^,a) = C(i,0)Ai(C(t,0),0). 



exp(~ I a-i£,{y')dy') = exp{a" / ^(C(s, a)) -ds) 







P 



hence the limit when a ^ of this quantity is 1 . 

For /3o given and ^0 = we identify ao such that, for a < and C < \i 
the asymptotic series defining Zj,mp converges to an analytic function. As 
C(t*,0) — if one chooses < ^ < jji, that is < < there exists 
ai > such that, for a < ai, ({t^,,a) < ^ < 

For t e ^^^d a < ai, C(i, a) < hence the functions (C(~7^))^-2p(— ^) 

and (^(^^))~^™(j(~7^) E^re well defined through the expansion of section 
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[31 Moreover, for t > ^ there exists C(^) such that (PT|) holds. Hence for 
t G we have the hmit of the functions Rp and Lk when a ^ 0. We are 

now ready to prove our main Theorem: 

Theorem 4 Let M he given. There exists a* > such that, for < a < a*, 
fi G [jj, M\ and \^\ < M, 5R7 > 0, the Evans function of the system Ev{a, f], 7) 
does not vanish. 

Proof Recall that we proved that 

Ev{a,(3,^) = EvoiP,j) + a-£'ui(/3, 7, a). 

The value of Evo{f3,j) (which is expressed through (|112p . (|114p ) docs not de- 
pend on t^,. Letting the leading order term of Evo{P,j) go to when t.^, — > +00 
(which gives, of course, the value of Evo{P,^) because it does not depend on 
t^) yields r = 1 as only possible positive solution for Ev{0, 7) =0 (see (|126p ). 
However, for r = 1 the remaining leading order term in t^ of Ev{0, P,"f) is not 
zero (see ([T^T])), hence a contradiction. 

The system of Kull-Anisimov has no bounded complex growth rate 7 
when Fr — 0{j) in the Hmit e ^ 0. To be more precise, recall that 

Evia,Pn) = e-^'+'^H^ eM-Io^ a^S,{y')dy')x 

[zi{Ri ~ jrLi - L2) + ^Z2{R2 + fjLi - U) 
Hz3(R3 + jL2 - jiU) 
+Cz4(i?4 - ^^3 - j^Lb - ie) 

+ rz|-("^iii + "^2-^2 + m^L^) + Y^{miLi + TO5L5 + meie)]- 

(112) 

We introduce Ap(C, a), Bq{(l, a), l<p<4, l<g<6 the hmctions given by 

/3ap(C, a) = MC. a),PCBg{C, a) = Bg{C, a). (113) 

We proved in Section[3]that the solution w^_^'' of (|42| which behaves as e^^^ ('^)-°'P)y 
when y +00, satisfying the condition is given through the relation (|59p 

^jfj~^J2p=i "^pdp — e~"^^"^^T(^)?i)^\ where {zj,mp) are given by Propo- 
sition d] through 



Zl = 


= /3 + /3ai(C,a) 


Z2 = 


= 7-/3 + /3H2(C,a 


Z3 = 




Z4 = 


= -/3 + /?|i4(C,a) 


TOl 


= /3fSi(C,«) 


TO2 


= /3p2(C,a) 


TO3 


= /3|S3(C,a) 


TO4 


= /3^B4(C,a) 


7715 


= /3^B5(C,a) 


7716 


= /3fS6(C,a). 



(114) 
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Equality ([55]) for a t such that ^(^1 Q^) < "^j along with f = a i ({t, a) i yields 

[(1 + CAi){Ri - j^Li - L2) + (e(r - 1) + C^)(^2 + ^ii - ^4) 
+ j^[BiLi + B2L2 + BsLs + BiL4 + B5L5 + BeLe]. 



We proved in Section [5] that the unique solution Uo{t, a) of the system (|85p 
satisfying the uniqueness condition ([90|l is 

4 6 

Uo{t, «) - ^ a) 4^ + 9g{t, 

p=l q=l 

and the functions i?j and Lp given by ((56|l satisfy the estimates 

4 6 

j=i p=i 

Note that the limit of the quantities i?" (as well as Rj when t +00) is known. 
We now consider the equality on Ev{a, (3, 7) when a ^ 0. The right hand side 
is independant of t because the left hand side is independant of t. Hence its 
value can be considered at t.^,. Once t, is fixed, we get the limit by taking a = 0, 
hence, using = — r — 1 



[(1 + CA\{0){R\{U) - LliU)) + CAl{Rl + rL? - L") ^15) 

+ai{Rl + rLl) + QAl{Rl-rLl~Ll) 

+C[5?L? + SOL° + Bill + BlLl + Bill + BlLlW 

where the relations are written at t = and at ^ = C(^*j 0) = 
Let Apj and Bq j being given through 

00 00 
KAl{0 = /3 5](/3CFAp^„/3CS°(C) = /3^(/3C)^B,,,. (116) 

By keeping only the leading order term in ^(y) for each equation in the system 
(which means that we consider the order of each quantity Zp and Mg), we 
obtain the recurrence system (jll7p : 
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vij + l)^ij+i = - rAi^j - 

+ 1) + l)A2j+i = rAaj - - - Sg^ 

+ 1) + l)A3j+i = Aij + + rAsj - Ai,j - B2,j - Be,j 

+ 1) + l)A4j + l = -Agj + rA4j + 

+ 1) + = -rSij - - Bij - rB^^j - + (r^ - 1)^2^ 

+ 1) + l)B2.j+i = Bij - rB2j + - rBe,j + rA^lJ + (r^ - 1)^3^ 
+ 1) + 1)53 j+i = S2,, - rSgj + Be J + (1 - r2).l4,, 
+ 1) + 2)B4j+i = ~Bi J + B5J + rA2,j + 
+ 1) + 2)B5.j+i = B4j - Bgj - A4J 
+ 1) + 2)B6.j+i = B3.J + Bs.j- + rA4j. 

(117) 

It is easy from this system to deduce that, under the hypothesis |r| < M , 
there exists a constant C > such that X]p=i l^p.il + — (t")''^!' 

hence ensuring that the analytic expansions defining A^ and i?° are extendible 
for all C,. The study of this recurrence system is the aim of the next paragraph. 



6.2 Behavior of the equivalent solution 

In what follows, we study the recurrence system. 
Consider , r) given by 

So(C, r)ii ^ -rii + ^i^ + 14 
Bo{^,r)i2 = is 
Bo{£,,r)i-i = ^zi - i2 + h 
Bo{£., r)ii = rii + (1 - r'^)ii - ji2 + 1*3 
Bo{^,r)i5 = is - C«4 + r£,ii. 

associated with the differential equation ~ BQ{S_,r)Y'^ . It is easy to check 

that B^^ is associated with the system (jllSp . obtained from ([61]l by taking into 
account the behavior of Zp and Mq that we obtained in Section [31 

dZ 



dt 
dMn 

dt 
dMfi 

dt 



We introduce 



^ =rZ2-Z3-(Afi+^M5) 

Z2 + - Z4 - {M2 + CMe) 
-Z3 + rZ4 + M3 

-r-Mi - M2 - e(M4 + rMs) - i)Zi + (r^ - 1)^2 
Ml - rAf2 + M3 - r^Me + (r^ - 1)^3 

i^=M2-rM:i+ ^Mq + (1 - r2)Z4 



dZ4 

dt 

dMi 

dt 
dM2 

dt 
dMj 



(118) 



-iMi 
M4 - Me 



M5 



.^3 



tM3 + Ms + ?Z4. 



i^o(2/) 



7C''-\h-9i+rg2). 



(119) 
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We have 

Lemma 7 Let H be the unique solution going to zero when y ^ of 

The coefficients Apj,Bq,j are the coefficients in the expansion in e(^(2/))~'^ of 
H, where ^{yy = vy. The vector (CA^, Ci?°) is equal to i?(ff ) for ( ^ ^(y)"". 

We need to obtain the relation with the asymptotic expansion in a(3 of 
Section [31 For this purpose, we introduce t = ey. We consider the following 
function, which will be tentatively a solution of (|118p 



(120) 



MM = j:,>iB'q^^eHay))-'''-\<i^ 1,2,3 

I Miiy) = Bl .e^{ay))-'''-\l = 4, 5, 6. 

Consider for example the first equation of (|12Qp . It rewrites 

hence 

^ J=0 j>l 

If we want to obtain {Apj,B'^j) indcpcndant on ^{y), it is a natural choice to 
write 

-# = r+^ (121) 

ay 

and, up to a change of origin in y, we obtain 

ay)-" = ^V- (122) 

Even with this equation, there will still be an additional term in the relation, 
related with j = 0. The resulting equation on Zi{y) given in (|120p is 

[Zi(y) - eA;,i(e(y))-n = C^a - r^Z^ - ^M,. 

ay 

With these two relations, we obtain the recurrence relation 

z.(j + l)A[^,+i = {A',^, - rA',^, - S^,), J > 1. 
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The same method applies to all the equations. 
If we impose the initial conditions: 

A' — _r+i A' — 4' — n 4' — 1 

nl 1 r)l r rtl rtl rtl pj _ n V^^'-'J 

^1,1 - --;;+iy^2,l - mi:^^3,l - ^4,1 - ^5,1 - ^6,1 " ^ 

which corresponds to the source term Fo{y) given by (jll9p . we have the identity 

^p,, =s; ,,.,Vj> l,Vp= 1,2, 3,4, V<z= 1,2,3,4,5,6 

where (A'p^^B'^j) are the coefficients of the expansion in e of the solution of 
(|118|) with the source term (|119p whereas {Apj, Bqj) are the coefficients of the 
expansion in C of {A^, B^). Lemma [7] is proven. 

Reduction of the problem Let us study the matrix Bq. Its eigenvalues are 
0, 1, and —1 and that associated eigenvectors are 

eo = (e,-2,0,0,l),ei - (1 - r, -i 0, 1, 0), Fi = (^,-1,1,0,1), 

e_i = (r+l,-i,0,-l,0),F_i = (^,-1,-1,0,1). 
The inverse matrix is given by 



' 2ii 


= ei ^ 


- e-i - 




-2eo) 




= Fi- 




-2eo 






= Fi - 










-(n 


- r)ei 


+ (r - l)e_i + 




. 2^5 




f F-i 


-e(ei +e_i). 





We rewrite a solution 1/ of 4^ = e-Bn^ as 

V = Foeo + Viei + WiFi + V-iC-i + W-iF_i. 

We obtain 

dV dVo dVi dWi ^ dV-i dW-i ^ „^ ,d^ . df_2,,, 

dy dy dy dy dy dy dy dy 

hence the associated system is 

' ^=e^[rHVo + W,+W.^)+rHVl+V^^)] 

< d^^eW,-el^[r\Vo + Wi+W-,)+r'{V, + V-,)] (124) 



49 



From the relation 



4(Ai/i + + + + Bigi + 3292 + B^g^ + B^g^ + B^g^ + Bege) 

= Ai{ei + e_i + ii?o) A (ei - e_i) + (rAz - Si)i?o A (ei + e_i) + A2R0 A (ei - 

+ {rA3 - B2){Fi ~ A (ei + e_i + ii?o) + ^3(^1 - A (ei ~ e_i)+ 

(rA4 + BsKei + e_i + ii?o) A (Fi - F_i - ^(ei + e„i)) 

+A4(ei + e_i) A (Fi + F_i - ^(ei + e_i)) 

+BiRo A (Fi - + B5R0 A (Fi + F_i - ^(ei + e-i)) 

+Be{Fi - F_i) A (Fi + F_i - ^(ei + e_i)). 

we deduce a new basis of A^(]R^) in which the coefficients are Ai, rA2 — -Bi, A2, 
rAs — B2, A3, rAi + B3, A4, B4, B5, Bq. The source term (|119p in the basis of 
A^(E^) associated with cq, e±i, F±i we find 

My) = -"^r^liiei + e_i) A (ei - e_i) - |(ei - e_i) A (Fi + F_i - 2eo)] 
+ (I;TT)^"''"'[(^i - ^-1) A (ei - e_i) - i(ei + e_i) A (Fi + F_i - 2eo 

(125) 



End of the proof The coefficient of ei A Fi in the source term is thus 

4 v + 1 

The theory of Fuchsian systems (see Hartmann [11]) shows that there ex- 
ists a constant a* such that the projection of H{y) on ei A Fi behaves as 
_ i^a.g2t^ri4 _ --i^)ei A Fi . The leading order term in of the Evans func- 
tions Fw(0,/3,7) writes 

Ev{0, f3, 7) = e('-+i)**iFC* 

- FO) + C(i*, 0)[i?(i?? - F") + i^(FO + rL? - L^) 
+AURl + rL^2)+Ai{Rl-rLl-Lt) 
+B1L° + B°,L° + B^.Ll + BlLl + F^ig + FgFg]] 

where iO(t*) = A%{t,,Q))e-^''t-''' , mu) = FO(C(i* , 0))e-2*-C"- • Hence, 



p 
as 

-ei A Fi = -^/i + f2 - h + h + rgi + (1 - r)(.g2 + 53) ~ ^(f4 + 55) 
and the hmit of (Fp, F°) when ^ -l-oo is given by ([SU)) . the Umit of 

F.(0,/3,7)(e('^+^)**iFC^^/3)-^ 

is 

(.-i)(-!:±l + -l^). (126) 
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As we seek positive values of r, the factor ei'>~+'^)t'- goes to infinity when — > 
+00, hence this limit is necessarily 0. The two possible values of r are thus 
r = 1 and r — hence r = 1. 

We also check that, for r — 1, the projection of Fo(?/) on the space associated 

(2) 

with the eigenvalue +1 of Bq ' , space generated by ei A eg and Fi A eg, is 

(t:^'^"^^ + 2(t/+i) ^~''~^)'^i ^ "^o- ^^'^'^ ^'^'^ relations 

ei A eo A S" 0, Fi A eo A S" = 

because eiAeo = -^/i+2/2+/4+gi- 1.95 andi^iAeo = -^31-^32+254+35+36- 
The leading order term of the projection on this space is 0. We get eiAe_i A5 = 
—Ail A Z2 A 13 A 14 A 15 because ei A e_i — — 2/i + |/2 + |3i, and similarily 
ei A F_i A 5" = 0, e_i A Fi A S" = -2ii A 12 A is A «4 A 15, Fi A F_i A S = 0. As 
the vector Fo{y) has the following projection on the eigenspace associated with 
the eigenvalue 

Pr{Fo{y)) =-2i;r'^(2eiAe_i + le_iAFi-ieiAF_i) 

+ ^(e_i A Fi + ie_i A Fi + (1 - i)ei A F_i) ^'^'^ 

the associated leading order term gives a non-zero contribution, hence a contra- 
diction. The main theorem is proven. 



7 Annex: Volterra type expansions of solutions 
of differential systems 

We consider a solution of associated with the growth rate A_(l) — a/3, 
which is the eigenvalue of smallest real part — Mq^-'(-I-oo). We recall that there 
exists a regular function //(a) such that 

A_(l)-a/3 = -l + a^(a). (128) 

We will prove that there exists a unique solution w^^ of ([33]) satisfying (|49l) . In 
order to prove the existence and uniqueness of w\_ we prove the following 

Proposition 6 Let U be the solution going to (1,0, ...,0) at +00 of the model 
system: 




= \,u, + (1 - ay)) E NAay))Uk{y), i<j<d 



where the properties of the complex numbers Xj and of the functions Nj^. are the 
following 
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= Ai <nX2...<^Xd-i, Ndd = 0. 
The function U is given by 



U{y)^ {1,0,.. ..,0) + {l-ay))w{y). 



Note that, in the hypothesis, N^d = is only there for simpUcity purposes 
and is obtained by considering in the last equation (in which we should have 
{l — £,)Ndd = ^'^^$t) the conjugation by the exponential of the primitive of 

Construction by recurrence of the Volterra operators We prove Propo- 
sition[5]by recurrence. Consider the last equation (line d of the previous system). 
We have 

^^x,Ud + Y.{i^ay))NdkUM- 

y fe=i 



This equation is equivalent to 

, d-l 

-{Ude'^-y) = J2''~^"'^^~^(y^^^dkU,{y). 



As Ud is bounded, the limit of UdC ^''^ is zero, otherwise Ud would not be 
bounded. Integrating from y to ya and letting yo go to infinity, we have 

d-i /•+00 

Ud{y) = / (^(^) - ^)Ndk{m)e-^^'U,{s)ds, 

k=i -^y 

and defining K^"^ through 

K^dfiU) - (1 ~ ay)r'e^'' r°^{^{-^) l)NdM^))e-^''U{s)ds (129) 
we obtain 

c/.(y) = (i-e(y))E<^(f^^)(y)- (130) 



fc=i 



Let us study the properties of k'^^ . 

Consider U such that, for ^ > (and yo such that — C(yo)) we have the 
estimate 

3N,C,yy>yo,\U{y)\<Cil~ay))''- (131) 
We obtain easily the estimate (thanks to \d-p > 0) 

yy > 2/0, \Ki]\u){y)\ < ^^^i^^(l - ^(2;))^. 
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Replacing (|130p in the rf — 1 first equations of the differential system of Propo- 
sition [6l we obtain 

1 < J < d - 1, ^ = A, [/, + (1 - ay)) E (^fc) 
^ fe=i 

where 

^^^(t^fe) = ^./c-C/fe + (132) 
We get the estimate, under the assumption (|13ip on U 

V > 2/0, \N^l\u)iy)\ < MoC'il + )(1 - ^(y))^. 

This rewrites 

J/>Z/o,|iV«(C/)(|/)|<AfiC(l-^(y))^ 

where Mi = A^^o(l + , Z'.'^" +i ). The procedure proceeds as follows: 
Let p be an element of 0..d — 2. We write the sequence of operators K^J^J'l 



d—p,k'! 

1 < fc < d - p - 1, and ivj^^^\ j < d - p ~ I, k < d - p - I such that 



d-p-l 

C/,_^(y) = (l-^(y)) E (133) 

^ = A,L/, + (1 - e(2/)) 'ivir'V.)(2/) (134) 
^ fc=i 



where 



^j^+i)([/,) = 7v(^)(L/fe) + (1 - ^(y))-^iv(l^((l - OKt'piiUk)) (135) 
and the operators K^f^^^, are constructed as follows using the auxiliary problem 

^ = A,„,F + (1 - C(y))^i!,),,,_,(l^) + /. (136) 
Bounded solutions when y —> +00 of (jl36p are obtained, for d ^ p > 2, through 

r+oc r+oo 

Viy)e-^--^y^ ias)-l)e-^'-'''N^/X^^_^{V)is)d.s- f{s)e-^--^^ds. 

■ly -ly 

This can be written, using g{y) = V{y)e~^'^-py , under the form 

g = Kd-pig)~ / f{s)e-^^--'ds. (137) 
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Hence we may write, for / satisfying (|13ip 

1=0 •'+°° 
which defines Td-p through 

V = Td-p{f). (138) 

Replacing / by (1-^(2/)) J2k=i ^d'-p.kiUk)iy), we obtain the expression of Ud-p 
in function ofC/fc,l<fc<d — 1 as 

Ud-piy) = Td^pi{i-o K-p.km) 

k<d—p 

hence 

(1 - avMt'piiUk) = r,_,((i - ONjfX.kiUk)). (139) 

Construction of the solution The last step of this recurrence is to construct 
the solution Ui. The equation on Ui writes then 

^ = (l-^(y))iVr^)(C/i), (140) 

where n['^~^'' satisfy the relations of Proposition [7l This is equivalent to 

Ui{y) - C/i(yo) = /'(I - m)N[t'\ms)ds. 

Using the limit Ui{yo) — s- 1 when j/o — > +oo, we obtain 

Uiiy) = l+ r (1 - m)N[t'\Ui){s)ds. 

J+oc 

We introduce the operator such that J7i = 1 + Ki{Ui). For g satisfying the 
estimate (|13ip . a consequence of proposition [7] is that 

\K\9)iy)\<^^,:^^^il-ayr'-'. 

Introduce Ui{y) — 1, and the sequence — K^{Ui^^). It is straightforward 
to show the inequality for y > yo 

\u^{y)\<{^r^,{i~ay)f 

hence as the series J2'n=q(7^)^ W^^ ~ Hv))^ normally convergent for y > 

so 

yo, the series X]^o iy) normally convergent and is the only solution going 
to 1 as y goes to +00 of (|140p . We have the estimate 

|[/i(y)|<exp(-^). 

so 

which is of the form p3T)) for = 0. 
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End of the construction As Ui satisfies an estimate of the form (|13ip for 
= 0, from the definition of the operators Kjf^^l we deduce that 



U2iy) = {1 - ay))w2iy) 

-y). Replacing in • 
inequality on [/i and on U2 wc deduce that 



where |w2(?/)| < ^^exp(-^). Replacing in the equality (|133p forp = d — 3 the 



M 

C/3(y)<M(l-e(y))exp(-^). 

so 

By recurrence on p we obtain the inequalities on Uj for j ^ 1, which proves 
Proposition [S] 



Estimates We prove in this Section the following 

Proposition 7 The o\ 

for g satisfying UlSl]) : 



Proposition 7 The operators N^^ and K^_^J'l satisfy the following estimates 



Vy > 2/0, \T{g){y)\ < -^^^—(1 - ^(y))^. 
This is a consequence of the more precise proposition 

Proposition 8 Under the assumption g satisfy the estimate I1131\) . the opera- 

(v) 

tors N^l, satisfy the estimate 



N 



\N\'^{g){y)\<CM,{l~av)) 
where Mp+i = Mp{\ + exp(-^^)) and 

so so 

We prove the second proposition by recurrence. The first estimate that we have 
to deduce from the recurrence assumption on n'^^ is the estimate on K^J^^\- 
To obtain this result, we have to study the behavior of K^-p through the esti- 
mate on N^^K We have, for / satisfying (|13ip for all N, the inequality 

/•OO 

\Kd-p{f)(y)\ < / {l^as)f^'Me^''-''^y-^Us. 
As iV + 1 > 1, we use y — s < and Xd-p > to obtain 
yy>ynAKd^My)\<M / [i - ^[s)r+^ds < -^^——{l ^ ^{y)r+\ 

■Jy So + ^) 
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Note that we increase the power of (1 —■?(?/)) in the result. 

The second estimate is based on the expression of T^-p obtained through (|138[) . 

For A'' > 1 and / satisfying (|13ip . we have 

The estimate on Kii_p, deduced from the recurrence hypothesis is 

-'y •'y so 

Hence we obtain 

hence the series X]i>o ^'^'^^^^ {^d-pY {J^ f{s)e~-^'^~p'^ds) is normaUy convergent. 
It defines the function Td-p{f) and we have 

\T,^p{f){y)\<C{\-ay)feM^)- 

so 

Finally, using this estimate on Td~^p as well as the relation (|139p , for / satisfying 
(|13ip we have 

1(1 - i{y))K'f^;i{f)\ < ai - e(y))^+^ exp(-l^) 

so 

which gives 

So 

Using the estimate on -^j^-p obtain 

1(1 - C)-^41p((l - < MpCil - ay)f exp(^) 

so 



hence 

\N^/,-^'\f){y)\<C{Mp + MpeM^m-ay)r 

so 

If we introduce Mp^i = Mp{l + exp(-^)), we thus deduce 

\N^l-''\f){y)\<CMp+,il-ay)f- 

The proposition [5] is proven. As we used the assumption that ^Xd-p > to 
ensure that a solution of (|136|) is given through 
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hence 

V{y)e-^'--y-V{y,)e-^---y' = f e-^^--'{l-as))N]l'}^^^_^{V){s)+S{s)e-^^--'ds. 

If Ve~^'^-'py has a hmit, then this is in contradiction with the fact that V is 
bounded for 5ftAd_p > 0, and this argument is no longer vahd if \d-p = 0. 
Hence this proves that the recurrence stops at p such that Xd-p — hence 

p = d — 1. The estimates of Proposition [7] are vahd for -/vj^^^\ hence the 
recurrence proceeds till SRA^-p > 0. This recurrence processus stops for \d-p = 
because we cannot assert that the equation (|136p has bounded solutions going 
to a constant for y +oo. 
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